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We show that off-site processes and multi-orbital physics have a crucial impact on the phase diagram of 
quantum gas mixtures in optical lattices. In particular, we discuss Bose-Fermi mixtures where the intra- and 
interspecies interactions induce competing density-induced hopping processes, the so-called bond-charge inter- 
actions. Furthermore, higher bands strongly influence tunneling and on-site interactions. We apply a multi- 
orbital interaction-induced dressing of the lowest band which leads to renormaUzed hopping processes. These 
corrections give rise to an extended Hubbard model with intrinsically occupation-dependent parameters. The 
resulting decrease of the tunneUng competes with a decrease of the total on-site interaction energy both affecting 
the critical lattice depth of the superfluid to Mott insulator transition. In contrast to the standard Bose-Fermi- 
Hubbard model, we predict a large shift of the transition to shallower lattice depths with increasing Bose-Fermi 
attraction. The applied theoretical model allows an accurate prediction of the modified tunneling ampHtudes 
and the critical lattice depth both recently observed experimentally. 
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Quantum gas mixtures in optical lattices are well suited to 
study in detail interaction-induced effects in condensed mat- 
ter. They allow for the investigation of systems with spin de- 
gree of freedom and with different species of particles that 
can even obey different quantum statistics. In particular, the 
experimental realization of atomic mixtures of bosonic and 
fermionic particles (e.g., ^^Rb-^°K) in optical lattices [1- 
4] triggered a vivid discussion on the role of inter- and in- 
traspecies interactions. These experiments allow for the ob- 
servation of the bosonic superfluid to Mott-insulator transi- 
tion in the presence of fermionic atoms. The prominent fea- 
ture observed in all experiments is the decay of visibility and 
condensate fraction of the bosonic subsystem induced by the 
interaction with the fermionic atoms. Two possible explana- 
tions for this drop in bosonic coherence were proposed. First, 
the process of adiabatic heating while ramping the lattice has 
been suggested [5, 6]. It is caused by different contributions 
of the atomic species to the total entropy and is therefore spe- 
cific to the loading procedure of experiments with ultracold 
gases. Second, an interaction-induced dressing of tunneling 
and interaction processes has been found that causes a shift of 
the superfluid to Mott-insulator phase transition [3, 7-9]. The 
latter effect corresponds to a necessary extension of the Hub- 
bard model at zero temperature and is therefore fundamental 
for various lattice systems. The important role of interaction- 
induced processes in optical lattices is caused by the specific 
shape of the Wannier functions and the possibility of high fill- 
ing factors. 

The standard Bose-Fermi-Hubbard model [10] is restricted 
to the lowest single-particle band and on-site interactions. In- 
terestingly, it fails to describe interaction effects in boson- 
fermion mixtures. Assuming a fermionic band-insulator, 
which is present in the experimental realizations [1-4], the 
boson-fermion interaction gives rise only to an irrelevant shift 
of the global chemical potential. Even for realistic assump- 
tions for the confining potential, the interspecies interaction 
has little influence [11]. Due to an effective screening even 
a small shift in the opposite direction is predicted. In con- 
trast to the experimental results, the superfluid phase is there- 
fore more stable within this framework, which poses the ques- 



tion of the applicability of the standard Hubbard model in this 
case. It was pointed out that the Bose-Fermi-Hubbard param- 
eters, i.e., tunneling and on-site interactions can be strongly 
influenced by the inclusion of higher orbitals and non-local 
interactions [7-9]. Off-site interactions include the so-called 
bond-charge interactions that have a direct density-dependent 
influence on the total tunneling [9, 12-18]. For purely bosonic 
systems and without taking off- site interactions into account, 
the inclusion of higher bands has been discussed in, e.g., 
Ref. [19-29]. Both extensions can be included in extended 
Hubbard models [9, 17, 18] and can be illustrated as effec- 
tive modifications of the tunneling potentials for both atomic 
species [17]. 

Here, we present quantitative results for the phase diagram 
of Bose-Fermi mixtures using an extended Hubbard model 
with multi-orbitally dressed processes [17, 18]. It includes 
density-induced tunneling, so-called bond-charge interactions 
[9, 12-18], as well as a multi-orbitally renormalized tunneling 
and interaction processes [9, 17, 18, 27]. The physical effects 
discussed here are in general present for all interacting quan- 
tum gas mixtures. The exact results, however, depend on the 
quantum statistics of the particles and the specific parameters 
such as the detuning from the light potential and the atomic 
masses. 

As a central result, we present the phase diagram of the 
superfluid to Mott-insulator transition in a Bose-Fermi mix- 
ture in section I. Furthermore, we discuss the crucial effect 
of off-site interactions and the corresponding shortcomings of 
the standard Bose-Fermi Hubbard model in section II. Subse- 
quently, we present the procedure to incorporate higher-band 
processes in section III. Afterwards, the corresponding ex- 
tended Hubbard models and the implications for the bosonic 
phase transition are discussed in detail in section IV. 



I. PHASE DIAGRAMS 

We will now first discuss the resulting phase diagrams of 
the bosonic superfluid to Mott-insulator transition in the pres- 
ence of a fermionic band-insulator, where the individual cor- 
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FIG. 1 : a Phase diagram for the superfluid to Mott-insulator transition of bosons with a fermionic band-insulator at different interspecies 
attractions. The predictions of the standard Hubbard models are shown as a dahed black line. The attractive interaction effectively reduces the 
total tunneling resulting in extended Mott- lobes, b The critical lattice depth of the superfluid to Mott-insulator transition as a function of the 
interspecies scattering length asF- The transition occurs at significantly shallower lattices than in the purely bosonic system (asF = 0). The 
dashed lines correspond to the Bose-Fermi Hubbard model and the dotted lines to an extended model with only lowest-band processes (section 
IV). 



rections to the standard Hubbard model are discussed in de- 
tail below. For concreteness, we choose a mixture of bosonic 
^''Rb and fermionic ^^K in an optical lattice with a spacing 
of a = 377 nm. For the respective wavelength, the Wannier 
functions of both species are almost identical. The interaction 
between the bosonic atoms is fixed to a repulsive scattering 
length of aBB = 102 ao [27], while the attractive interaction 
between the two species is tunable over a wide range using 
a Feshbach resonance [3]. The fermionic nature of the spin- 
polarized potassium atoms simplifies the system, as we can 
assume a band-insulating phase and thus a fixed atom-number 
of one fermion per lattice site. Note that it is possible to di- 
rectly apply the presented methods and extensions to a variety 
of atomic quantum gas mixtures in lattices. In particular, both 
atomic species can be bosonic or fermionic and the generaliza- 
tion to multi-component systems with more than two species 
is straight forward. 

For a fermionic band insulator, the fermionic degrees of 
freedom are frozen out and the physics can be described by 
an effective bosonic model that takes into account all effects 
induced by the interaction with the fermions. In the frame- 
work of this paper, we will discuss in detail the derivation of 
an effective Hamiltonian which reads 

^ext = - E ^jJ^^n. + ^ ^n. - /i E 

We will see that despite its simplicity it already includes 
higher-band and bond-charge off- site processes. The latter 
gives rise to an occupation-dependent tunneling = 
-^B + (^i + — 1)Xbb + 2Xbf even within the lowest 
single-particle band. Here, Jb is the conventional tunneling; 
Xbb and Xbf are the bond-charge tunneling elements aris- 
ing from Bose-Bose and Bose-Fermi interactions, respectively 
(section II). The interaction induced occupation of higher or- 
bitals leads to a further occupation dependency of all parame- 
ters, i.e., Jb, Xbb, Xbf and En. 



The tilde above the parameters and operators in (1) indi- 
cates the multi-orbital dressing as discussed in section III. The 
effective single-band Hamiltonian (1) uses the ground state of 
the interacting system, called the dressed band, instead of the 
lowest single-particle band [17, 18]. The dressed operators bi 
and b] annihilate and create bosonic particles on site i in this 
dressed band and hi = P-bi counts the number of bosons on 
site i. It is important to note that after the transformation to 
the dressed band, the phase diagrams can be calculated using 
standard single-band methods. The renormalized on-site en- 
ergy En is composed of the single particle energies of bosons 
eB,n and fermions eF,n as well as the interaction energies for 
the repulsion between the bosons ^n{n — l)Un and the attrac- 
tion between the species nUBF,n- The chemical potential fi 
fixes the total number of bosonic atoms. 

After calculating the dressed parameters, we apply 
Gutzwiller mean-field theory to compute the critical lat- 
tice depth of the transition from the superfluid to the Mott- 
insulator. The phase diagrams of the extended model (1) are 
shown in Fig. la. The effective chemical potential ji — Ei 
is given in units of the Hubbard on-site interaction U, where 
El is the renormalized on-site energy of one boson and one 
fermion. For vanishing interaction between the bosons and 
fermions aBF = and a repulsive interaction aBB = 102 ao 
among the bosons, the Mott-lobes are contracted compared 
with the standard Hubbard model. This is a result of a de- 
crease of the on-site energy and an increase of the total tun- 
neling caused by off-site interactions [17]. For increasing at- 
traction between bosons and fermions the effect is reversed 
[9] and the Mott-lobes are extended exhibiting a critical tran- 
sition point at much lower lattice depths. This effect can be 
attributed to a strong reduction of the total tunneling ampli- 
tude induced by off- site interactions. 

In Fig. lb the critical lattice depths for the superfluid to 
Mott insulator transition for one, two and three bosons per 
lattice site are shown as a function of the interspecies inter- 
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action strength. The soHd Hnes depict the results obtained 
using the extended model (1), while the dashed lines corre- 
spond to the standard Bose-Fermi Hubbard model, which pre- 
dicts no dependency on the interspecies interaction aer (sec- 
tion II). As discussed above, for gbf = the Mott-insulator 
transition is shifted to deeper lattices, where the shift is in- 
creased with the bosonic filling. This is mainly caused by the 
bosonic bond-charge interaction Xbb enhancing the total tun- 
neling. With increasing attractive interaction asF the transi- 
tion is strongly shifted to shallower lattice depths. Depending 
on the filling, the shift of the Mott-insulator transition caused 
by the fermionic atoms is 3-4 for qbf = —300 ao. 

The predicted shift is considerably larger than calculated 
with the adiabatic band elimination method in Ref . [9] which 
also incorporates the bond-charge interactions. The effective 
potential approach in Refs. [3, 7] does not include the im- 
portant contributions of fermionic on-site energy and bosonic 
bond-charge interaction. The extended Bose-Fermi-Hubbard 
model (1) discussed here contains all relevant energies that 
can affect the superfluid to Mott-insulator transition at zero 
temperature. As a result, the Mott-insulator shift in Ref. [3] 
can be partly explained by interaction-induced effects. This 
provides a consistent picture, where the experimental obser- 
vations [1-3] are a combined effect of the Hubbard extensions 
and the adiabatic heating processes [5, 6] which depend on the 
initial temperature of the quantum gas. 

II. OFF-SITE INTERACTIONS 

We now turn back to the full description of the applied ex- 
tended Hubbard model which features two corrections to the 
standard Hubbard model. First, off-site interactions lead to 
a significant contribution to the total tunneling amplitude by 
changing the effective tunneling potential. Second, the in- 
clusion of multi-band processes causes a modification of all 
model parameters. 

As mentioned above, mixtures of ultracold spin-polarized 
bosonic and fermionic atoms in optical lattices are usually 
described by the standard Bose-Fermi Hubbard model. The 
underlying tight-binding approximation restricts the model to 
the lowest single-particle orbital and interactions between par- 
ticles on the same lattice site. The resulting Hubbard Hamil- 
tonian reads 

^BFH = - E(^B6l6,+ Jf///,) + ^^Mn^ " 1) 

+ ^ UBFTlirhi - ^{llBTli + flFm), 
i i 

(2) 

Here, bi (fi) is the bosonic (fermionic) annihilation opera- 
tor and hi (m^) the respective particle number operator. In 
general, the tunneling matrix elements for bosons (Jb) and 
fermions (Jf) can have different values. The on-site inter- 
action is fully described by the parameters Ubb and Ubf 
for intra- and interspecies interaction, respectively. The to- 
tal number of bosonic and fermionic atoms are fixed by the 



chemical potentials /iB and /iF- Under common experimental 
conditions [1-4], the fermions are in a band-insulator phase 
where Pauli-blocking prohibits tunneling. This freezes out the 
fermionic degrees of freedom and the resulting Hamiltonian 
captures the behavior of the bosons under the influence of ex- 
actly one fermion per lattice site. Consequently, we can set 
fifj ^ 0, ^ 1 and get 

^FBI = -Yl -^B^l^i + ^ ^^(^^ ~ ^) 

+ X^(^BF - Mb)^*- 

i 

The interaction energy Ubf between bosons and fermions 
can be absorbed into an effective chemical potential /ieff = 
l^B — Ubf and the resulting Hamiltonian does not differ from 
the standard Bose-Hubbard model. Thus, the behavior of the 
bosons is not influenced by the homogeneously distributed 
fermions, which is in contradiction to the experimental ob- 
servations [1-3]. 

In the derivation of the standard Hubbard model, it is ar- 
gued that interaction processes between particles on neighbor- 
ing lattice sites can be neglected due to their small amplitudes 
compared with on-site interactions. This argument is how- 
ever only partly correct, since some of these processes, such 
as bond-charge interactions, involve the hopping of particles. 
Compared to the conventional tunneling, these processes can 
be non-negligible and consequently alter the phase diagrams. 
In particular, the Wannier functions in optical lattices differ 
strongly from their counterpart in solid state materials lead- 
ing to comparably large matrix elements for bond-charge pro- 
cesses. In addition, the possibility of larger fillings in bosonic 
systems can enlarge these interaction effects. Consequently, 
off-site interaction processes can be strongly enhanced for op- 
tical lattice systems. 

Consider the interacting part of the full two-particle Hamil- 
tonian for the lowest band and two neighboring lattice sites L 
and R 

Ant = \Y. Uf^i hmhk + J2 Ul^ki bUjfkbi, (4) 

ijkl ijkl 

with j, /c, / = L, R and 

,-xBB/BF f B*/ \ B/F*/ /\ t^/ /\ 

Uijki =9bb/bf < (r)^/ (r)y(r,r) 

J (5) 

xwf/^(r')wP(r) d^rd^r'. 

B /F 

The basis functions w- ' (r) are the maximally localized 
Wannier functions that describe a boson/fermion sitting on 
site i and the interaction strengths are given by ^bb = ^|~^b 

and ^BF = ^^^^^BF with the mass of the bosonic atoms niB 
and the reduced mass mr of boson and fermion. The inter- 
action potential ^ (r, r') we applied describes the scattering 
properties using a finite-ranged box potential (see [17]). In a 
lowest-band treatment, this is usually replaced by contact in- 
teractions, i.e., a ^-pseudopotential. In a multi-orbital frame- 
work the latter would lead to mathematical subtleties [25]. 
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FIG. 2: Interaction-induced off-site processes. Processes induced by the Bose-Bose interaction are a bosonic on-site interaction, b density- 
density interaction, c pair- tunneling, and d bond-charge tunnehng. Analogous interaction processes are induced by Bose-Fermi interaction, 
namely, e on-site interaction, f density-density-interaction, g cross-tunneling, h pair-tunneling, i bond-charge tunneling of a boson and k bond- 
charge tunneling of a fermion. The shading indicates the most important processes that are taken into account in the extended Hubbard model 
(1). The bond-charge tunneling processes (d and i) are illustrated in 1 and m via tunneling in effective potentials. 



The distinct processes arising from the full two-body 
Hamiltonian (4) for Bose-Bose and the Bose-Fermi interac- 
tion are depicted in Fig. 2. In addition to the on-site interac- 
tion, the Bose-Bose interaction leads to the density-density 
interaction process Vbb '^i'^j, the correlated tunneling of a 
particle pair Pbb ^l^^j. and the bond-charge assisted tunnel- 
ing — Xbb bl{ni-\- nj)bj. The respective matrix elements are 

^BB U^fi, Pbb ^^jV^' andXBB -U^Fj ~^^fj' 
The amplitudes of these processes are plotted in Fig. 3 as 
solid lines. While the density-density interaction Vbb and 
the pair tunneling amplitude Pbb are several orders of mag- 
nitudes smaller than the the standard Hubbard processes, the 
bond-charge interaction Xbb is only one order of magnitude 
smaller than the conventional tunneling Jb. As the bond- 
charge interaction — Xbb + ^j)bj scales with the par- 
ticle number on both involved sites, it can easily reach non- 
negligible values and must be accounted for (see also [17]). 
In addition, it scales linearly with the interaction strength 
which can be tuned experimentally. Here, it is essential that 
the bond-charge interaction contributes to the tunneling of the 
particles rather than the on-site interaction. 

For the Bose-Fermi interaction two more processes exist as 
the interacting particles are distinguishable (Fig. 2b). First, 
the cross-tunneling, which is similar to the density-density in- 
teraction except that the particles interchange, and second, the 
bond-charge interaction, where either a bosonic or a fermionic 
particle tunnels. However, assuming a fermionic band insula- 
tor all processes involving the hopping of a fermion are for- 
bidden. The respective amplitudes of the processes induced 
by Bose-Fermi interaction are plotted as dashed lines in Fig. 3. 

In conclusion, only on-site interactions and bond-charge 
tunneling of bosons have to be taken into account (shaded in 
Fig. 2), since other processes are prohibited by the fermionic 
band-insulator or contribute only with extremely small ampli- 
tudes. Thus, the necessary extensions of the Hubbard Hamil- 
tonian, i.e. the bond-charge processes in Fig. 2d and i, read 



X = X 



BB 



X 



BF 



fi fi ~^ ^ijjjfjfj^^j' 



(6) 



An intuitive physical understanding of the bond-charge tun- 
neling induced by Bose-Bose as well as by Bose-Fermi inter- 
action can be obtained by the analogy to an effective tunneling 
potential [17]. Assuming contact interactions in the single- 
band description and using the integral expressions (5), we 
can rewrite the expression (6) as 



X = h]L 



I' 



'(r) {gBBp- 



BB 



gBFpf)wf{r), (7) 



where we introduced the reduced densities 



ni\w. 



(r)|^ + (n, 
«;f(r)p 



l)kf(r)| 



B/ 

3 



(8) 
(9) 



The — 1 in the bosonic density corresponds to the exclusion of 
self-interactions and is directly obtained from the commuta- 
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FIG. 3: Matrix elements for the different hopping and interaction 
processes. The solid lines are the purely bosonic processes, while 
the dashed lines correspond to Bose-Fermi interaction and fermionic 
tunneling. Shown are the intra- and interspecies on-site interaction 
Ubb and Ubf (green), the conventional tunneling amplitudes Jb and 
Jf (red), the bosonic bond-charge tunneling Xbb and Xbf (blue) 
as well as density-density interactions Vbb and Vbf (orange) and 
correlated pair- tunneling Pbb (cyan). 
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tion relations. Inside the integral, we can use the locality of 
the Wannier functions to replace these operators by the density 
functions pf^ pb{y) -\wf{Y)\^ and pff py{v). The 
bond-charge tunneling operator (7) can now easily be unified 
with the conventional tunneling to find the expression 



J + X={wf{v)\ 



PI 

2m 



+ Kff(r))K(r))6j6,- 



(10) 



which corresponds to the conventional tunneling in an effec- 
tive potential yeff(r) = V{y) + ^BB(pB(r) - (r)P) + 
^BFPF(r). 

In Fig. 2k and 1 the tunneling in effective potentials is 
sketched. Repulsive interactions, as between the bosons, ef- 
fectively reduce the lattice depth, while attractive interactions 
have the opposite effect. Depending on the relative scatter- 
ing lengths, the total tunneling can be strongly enhanced or 
reduced. This is consistent with the results for the fermionic 
tunneling obtained in [4]. The modification of the total tun- 
neling in the lowest band already leads to a considerable de- 
formation of the well-known phase diagram of the superfluid 
to Mott-insulator transition for bosons as discussed later in 
detail. 



III. MULTI-ORBITAL RENORMALIZATION 

Whereas in the previous section we introduced off-site in- 
teractions as an important extension to the standard Hub- 
bard model, we will now discuss another important feature of 
the Hamiltonian (1), namely the effective inclusion of higher 
bands. In the standard Hubbard model approach, only the 
lowest single-particle band is assumed to be occupied. How- 
ever, in the strongly correlated system particles are promoted 
to higher orbitals due to the interaction induced coupling be- 
tween the bands. By changing their wave functions the par- 
ticles minimize their on-site interaction energy (see section 
III.B). In particular, this also affects the tunneling and other 
off- site processes as outlined now in section III. A. 



A. Multi-orbital dressing 

For the following calculation of the multi-orbital effects, 
we assume the orbital occupation to be determined purely 
by on-site interactions, while off-site processes can be ne- 
glected due to the much smaller amplitudes. The prob- 
lem can thus first be reduced to a single lattice site. The 
single-site ground state is a superposition of many-particle 
states ^(n) = Eiv,M c^,m |M). Here, |A^) |M) = 
|no,ni, ...) |mo,mi, ...) is the product state with bosons 
and TTio, fermions in the Wannier orbital w^^^{y), where a 
indicates the orbital. The state ^(n) consists of the low- 
est single-particle band dressed with small contributions of 
higher bands. At zero temperature the particles will exclu- 
sively occupy the dressed band ^(n) instead of the lowest 
single-particle band. 



In general, any multi-orbital two-site operator acting on 
sites L and R can be decomposed in operators of the form 



O 



(11) 



where the summation is over all possible sets of orbitals 
{a} = {^1,^2,...} and {/3} = {/3i, ^2, ...}• A^'W is the 
amplitude for the corresponding process and oj^^ consists of 

creation and annihilation operators b[^'^'^^ and ^-"''^ for par- 
ticles on site i in the orbital ak. In the simplest case - the 
conventional tunneling from the right to the left site - we have 
only sets with a single orbital {a} = a and {(3} = /3, the op- 



erators on the left and the right site Og = b[^'^\ = b\ 



m 

R ' 



and tunneling amplitudes A"'^ as defined below. Instead of 
collapsing this operator to the lowest single-particle band, we 
reduce it to the dressed band [17, 18] given by the many- 
particle ground state ^(n). The effective operator O in the 
dressed band takes the form 



= A 



(12) 



where Oi is obtained from Oi by replacing all creation and 
annihilation operators S^^''^^ and S^^''^ by their counterparts of 
the dressed band b] and bi. The operators of the dressed band 
fulfill the usual relations bi |^(n))- = v^|^(n — 1))^ and 



The effective amplitude 



is obtained from the ma- 



trix element (^f| O where ^i(nL, tir) denotes the ini- 
tial and = ^(^l',^rO the final state of the process. 
It thereby includes the summation over all multi-orbital pro- 
cesses. Since the states are product states of the individual 
lattice sites [^(til)) [^(n^)), also the effective amplitude A 
decomposes into individual site contributions 



1 

N 



J2 A^^^^^^U^M\dt^\^{nj^)) 



(13) 



xWV)loi^M^(nR)), 



where TV = (^f| OlOr is needed for the correct nor- 
malization. Note that the effective amplitude is intrinsically 
occupation-dependent. 

As an example, for the conventional single-particle tun- 
neling of bosons it follows A^^^ = —JaSa^(3 and = 



(«)\ 



IS 



V^V^^^TT. Here, Ja = - {w^'^^l ^ + V{r) \w 
the tunneling amplitude in band a. As another example, the 
multi-orbital bosonic bond-charge operator 



E 



(14) 



decomposes in left and right part oj^^^ and that consist 

of either one or three creation/annihilation operators. In con- 
trast to the conventional tunneling, orbital-changing processes 
are allowed for the multi-orbital bond-charge operator. 
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FIG. 4: a Contributions to the total tunneling in the lowest band 
(dashed) and the dressed band (solid) for one bosons and one fermion 
on each site. The total tunneling (black) is composed of conven- 
tional tunneling J (green), bosonic bond-charge tunneling Xbb 
(red), and fermionic bond-charge tunneling Xbf (blue). The ab- 
solute value of all matrix elements are enhanced by orbital degrees 
of freedom but due to their opposite signs the multi-orbital effects 
partly compensate each other, b Deviations of the total tunneling 
-^n?n = -^B + (riL + TiR - 1)Xbb + 2Xbf from the standard Hub- 
bard tunneling Jb for different bosonic occupations riL = riR. 



The renormalized amplitudes are depicted in Fig. 4 for 
all relevant contributions to the total tunneling in the case 
of one boson and one fermion per site (solid lines). The 
amplitudes for conventional tunneling JnL,nR (green line) as 
well as bond-charge induced tunneling XBB,nL,nR (red) and 
-^BF,nL,nR (bluc) Can differ strongly from the lowest band 
values, which are indicated as dashed lines. The multi-orbital 
renormalization can enhance the conventional tunneling by up 
to 30% and the bond-charge induced processes can even be 
twice as strong as in the lowest single-particle band. 



B. On-site problem and on-site interaction 

This section is dedicated to the explicit solution of the 
many -body problem on a single lattice site for a given number 
of particles. For the multi-orbital dressing procedure (section 
III.A) we assume the ground state of the interacting system 



to be a product state of these solutions. We apply the method 
of exact diagonalization to compute the ground state ^(n) of 
n bosons and one fermion on a single lattice site and thereby 
obtain the solution of the full lattice problem. In particular, 
this leads directly to an occupation-number-dependent on-site 
energy En- The respective single-site Hamiltonian reads 

OL Oi 

+\ E ^bb"'''^ u-mm-'W^ ^^^^ 

+ E C/^f'''^S(")t/(/3)t/(7)5W^ 

with the bosonic particle number operator n^"^ = 
and the single-particle energies The operators rri^^\ 

j(a)t^ f{oi) ^j^g energy ep^^ are the fermionic analogues. 
The multi-band interaction amplitudes are defined as 



(a/375) 

BB/BF - ^BB/BF 



W 



(r)<%*(rOnr,rO 



(16) 



The interaction potential F(r, r') is a finite-range box poten- 
tial with a width of 7.5 nm (see Ref. [17] for details). We 
expand this Hamiltonian in the basis of many-particle Fock- 
states \N) \M). For n bosons and 1 fermion, the ground state 
is found by applying an exact diagonalization, the so-called 
configuration interaction method. We restrict the calculation 
to the lowest 9 bands per spatial direction and use a high en- 
ergy cutoff [31]. The on-site energy is directly obtained as the 
lowest eigenvalue of the matrix and its contributions can be 
computed as expectation values of the individual operators in 
(15) using the corresponding eigenvector. 

The energy contributions are plotted in Fig. 5a for Vb = 
15 Eyi and n = 3 bosons (solid lines). The values significantly 
deviate from the lowest-band approximation (dashed lines). 
The single-particle energies nee + ep (green line) are mea- 
sured relative to the lowest-band values and thus are always 
positive. The occupation of higher orbitals causes a contrac- 
tion of the wave functions which leads to an increase of the 
absolute value of the Bose-Fermi interaction. This results in a 
large reduction of the total on-site energy (black line) for large 
scattering lengths asF- Of course, also the repulsive interac- 
tion among the bosonic particles (red line) contributes but is 
less drastically influenced. 

The total on-site energies for various boson numbers are 
shown in Fig. 5b. The dashed lines correspond to a calcula- 
tion with a three times larger interaction range. In general, 
the contributions of higher bands are reduced with an increas- 
ing interaction range. The figure shows that the energy is 
only weakly affected by a the change of the (finite) interaction 
range. Additionally, we applied scaling theory to estimate the 
value of convergence for the on-site energy at the most de- 
manding parameters Vq = 20 £^r and aBF = —300 ao for 
an interaction range of 7.5 nm. We see that the error (dif- 
ferences in energies) converges exponentially with both the 
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length of the many-particle basis and the number of orbitals. 
Scaling our results according to the exponential behavior, we 
are able to determine the converged energy value, where we 
first perform the scaling in respect to the basis length for a 
given number of orbitals. The results indicated as crosses in 
Fig. 5b show only small deviations and justify the constraints 
for the basis length and the number of orbitals applied for the 
solid lines. 

The total on-site energy becomes intrinsically occupation- 
dependent (beyond the dependency in the standard Bose- 
Fermi Hubbard model (2)) and can be written in terms of ef- 
fective n-particle collisions [18, 26, 27] 

En = nEi + in(n - 1)^2 + ^n(n - l)(n - 2)Es... (17) 

The first term Ei describes the interaction energy by two- 
particle collisions between bosons and fermions. The sec- 
ond term is the interaction energy caused by processes that 
involve two bosons E2 = E2 — 2Ei and the third term 
Es = Es — 3E2 — 3 El involves three bosons. Although, the 
restriction to the first three terms is enough to describe well 
the energies for up to n = 5 bosons, we use here the exact 
values for E^ . 



IV. SUPERFLUID TO MOTT-INSULATOR TRANSITION 




-200 
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After discussing off-site interactions and multi-orbital 
renormalizations, we now turn back to the full many-body 
quantum gas problem. Obviously, from the above results, it 
is necessary to take both bond-charge interactions and higher 
bands into account. First, we now define an extended model of 
the lowest band that includes off-site interactions and discuss 
its implications. Afterwards, we replace the lowest single- 
particle band and parameters with the respective dressed ana- 
logues and thereby include higher bands in a very efficient 
way. 

The extended Hubbard model of the lowest band reads 

^ext = -Yl [JB^XBBifli + flj - l)+2XBF]^l^j 

Ubb ^^^^ 

i i 

While the repulsive interaction between the bosons increases 
the total tunneling, the attractive fermions reduce the bosonic 
mobility. As one central result and in strong contrast to the 
predictions of the standard Hubbard model (see section II) the 
superfluid to Mott-insulator transition is shifted. The phase di- 
agrams are shown in Fig. 6 for different attractive Bose-Fermi 
interaction strengths. For strong Bose-Fermi attraction and 
low bosonic filling, the transition occurs at much shallower 
lattices, due to the effectively deepened tunneling-potential. 
The effect is reversed when the repulsion between the bosons 
becomes stronger than the attraction to the fermions, which 
is the case for weaker Bose-Fermi interaction and higher 
bosonic filling. In Fig. lb the critical lattice depth for the 



FIG. 5: a Total on-site energy En and its contributions for n = 3 
bosons and m = 1 fermion as functions of the interspecies scattering 
length. The contributions are single-particle energies eB,F (green), 
Bose-Bose interaction Ubb (red), and Bose-Fermi interaction Ubf 
(blue). The total on-site energy is visibly lower than the lowest band 
prediction (dashed lines), b The total on-site energies En for various 
numbers of bosons n and one fermion. With increasing interspecies 
interaction the total energy decreases non-linearly. The dashed lines 
correspond to a three times larger interaction range. The markers are 
obtained by applying scaling theory as described in the text. 

transition is plotted as a function of the interspecies interac- 
tion strength aer- The dotted lines correspond to the extended 
Hamiltonian restricted to the lowest band. 

When including higher bands we must replace the lowest- 
band operators with those of the dressed band bi and b]. Also, 
the parameters J, X and U must be renormalized as dis- 
cussed in section III. All tunneling contributions, i.e., con- 
ventional tunneling, both bond-charge interactions as well as 
their multi-orbital renormalizations can be combined to one 
total tunneling parameter 

JnT,nj = Jni,nj-^XBB,ni,nj{'rii-^nj—l)-^2XBF,ni,nj, (19) 

which is explicitly occupation-number dependent. The renor- 
malized on-site energy 

En = nCB^n + eF,n + ^^(^ - ^)Un + riUBF,n (20) 

is composed of the renormalized single particle energies of 
bosons eB,n and fermions eF,n as well as the interaction ener- 
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FIG. 6: Phase diagrams of the lowest-band model (18) with bond- 
charge interactions for different scattering lengths. For comparison 
the results of the standard Hubbard model are shown as a dashed 
black line. 

gies for the repulsion between the bosons ^n{n — l)Un and 
the attraction between the species nU^Y,n- This allows to de- 
fine the extended Hamiltonian of the dressed band in equation 
(1), namely, 

^ext = - E ^jJ^^n. + E - ^ E ^21) 

The dressed-band Hamiltonian now takes into account all 
higher-band processes and all relevant nearest-neighbor in- 
teractions. The multi-orbital corrections of the Bose-Fermi 
interaction have a strong impact on the chemical potential at 
which the transition to a certain Mott-lobe occurs. This dis- 
torts and shifts the phase diagram along the axis of the chem- 
ical potential. Since this effect is not of any physical interest, 
we plot Fig. la in terms of an effective chemical potential 
/^eff = l-i — El. Concerning the transition point of the bosonic 
superfluid to Mott-insulator transition, the reduction of the to- 
tal on-site energy by multi-orbital processes counters the ef- 
fect of reduced total tunneling. Nonetheless, the total effect on 
the transition can be a shift of several recoil energies de- 
pending on interaction strengths and filling factors (Fig. lb). 
Note that the lowest-band model H^^t underestimates the im- 
pact on the superfluid to Mott-insulator transition by only up 
to l^R, which is surprising keeping the strong changes of the 
individual amplitudes in mind. However, this (coincidental) 
compensation of the contributing amplitudes depends on the 
choice of system parameters. Furthermore, it has been demon- 
strated that the on-site interaction [27, 30] and the effective 
tunneling matrix element [4] are experimentally well accessi- 
ble and can be measured independently. 

V. DENSITY-DENSITY INTERACTIONS AND 
PAIR-TUNNELING 

In the above model, several off-site processes have been ne- 
glected due to their small amplitude in the lowest band. In this 
context the question arises, whether the multi-orbital dress- 



ing can enhance them to non-negligible values. The corre- 
lated pair tunneling and the density-density interactions both 
have very small amplitudes because in the integrand the small 
tail of the Wannier function enters quadratically on both lat- 
tice sites. By contrast, in the bond-charge integral the tail is 
multiplied three times with the center of a Wannier function. 
However, when taking into account strongly delocalized wave 
functions of higher bands, this argument no longer applies 
and all processes need to be reconsidered, since the overlap 
integrals become comparable for all types of off-site interac- 
tions. These contributions are strongly suppressed in the case 
of density-density interactions, because the initial and final 
state both depend on higher-band contributions with small co- 
efficients Cat <C 1 to produce a large integral, whereas in the 
case of pair tunneling one of them can be the ground state 
(cat ~ 1). Although, the multi-orbitally renormalized pair 
tunneling is usually smaller than the conventional tunneling 
it can, in general, reach the same order of magnitude. Due 
to this structure it has very bad convergence properties in re- 
spect to the total number of contributing bands. As a fourth 
order contribution in perturbation theory it can, however, be 
neglected even for rather large amplitudes. Density-density 
interactions are not as strongly influenced and remain small. 
This justifies, the restriction to the extensions incorporated in 
Hamiltonian (1), which includes all relevant off-site interac- 
tions. 



VI. CONCLUSIONS 

We have discussed the important role of interaction effects 
in atomic quantum gas mixtures in optical lattices. Off-site 
interactions as well as higher band processes turn out to have 
a strong impact on these systems, which we are able to calcu- 
late using an extended occupation-dependent Hubbard model. 
In particular, we have focused on Bose-Fermi mixtures in this 
paper, where the standard Bose-Fermi-Hubbard model fails 
to cover all relevant processes. This manifests itself in a 
strong shift of the superfluid to Mott-insulator transition in 
the bosonic subsystem, which is not predicted by the stan- 
dard Hubbard model. The critical lattice depth is shifted to- 
wards shallower lattices with increasing Bose-Fermi attrac- 
tion. Similar corrections are present for all experiments with 
optical lattices and can be expected to be relevant, e.g., for 
Bose-Bose mixtures, low-dimensional systems or other lat- 
tice geometries. Omitting the condition of having a fermionic 
band insulator, which we have applied here, the presented ex- 
tensions of the Bose-Fermi Hubbard model can lead to very 
rich physics such as the formation of polarons. 

We have shown that for optical lattice systems the bond- 
charge tunneling [9, 12-18] is the most important contribu- 
tion of the nearest-neighbor off-site interaction as it can dras- 
tically influence the tunneling. Repulsive interactions en- 
hance the total tunneling, whereas the attractive interactions 
reduce it. In an intuitive picture, this can be described as low- 
ered and increased effective potentials, respectively. Further- 
more, higher band processes not only reduce the total on-site 
energy but also have an impact on the conventional tunnel- 
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ing and the bond-charge interactions. We have treated the 
problem by dressing [17, 18] the lowest single-particle band 
with interaction-induced occupations of higher-orbital states. 
This leads to a renormalization of interactions and tunneling 
parameters that become intrinsically occupation-dependent. 
These parameters have been used in order to define an ex- 
tended Bose-Fermi-Hubbard model capable of describing ef- 
fects of higher orbitals and off-site interactions in Bose-Fermi 
mixtures. 

The results show in general that interactions in multi- 
component systems can have a crucial impact beyond the stan- 
dard Hubbard treatment. In the presented case, the standard 
Hubbard model is incapable of describing the interspecies 



interaction between bosonic and fermionic atoms correctly. 
While here mainly the effects on the bosonic atoms have been 
discussed, the mutual interaction affects the fermionic atoms 
similarly, which has recently been observed in experiment [4]. 



VII. ACKNOWLEDGMENTS 

We thank U. Bissbort and W. Hofstetter for stimulating dis- 
cussions and acknowledge financial support by the German 
science foundation DFG under grant FOR801. 



[1] K. Gunter, T. Stoferle, H. Moritz, M. Kohl, and T. Esslinger, 

Phys. Rev. Lett. 96, 180402 (2006). 
[2] S. Ospelkaus, C. Ospelkaus, O. Wille, M. Succo, R Ernst, 

K. Sengstock, and K. Bongs, Phys. Rev. Lett. 96, 180403 

(2006). 

[3] T. Best, S. Will, U. Schneider, L. HackermuUer, D. van Oosten, 
L Bloch, and D.-S. Luhmann, Phys. Rev. Lett. 102, 030408 
(2009). 

[4] J. Heinze, S. Gotze, J. S. Krauser, B. Hundt, N. Flaschner, D.-S. 
Luhmann, C. Becker, and K. Sengstock, Phys. Rev. Lett. 107, 
135303 (2011). 

[5] M. Cramer, S. Ospelkaus, C. Ospelkaus, K. Bongs, K. Seng- 
stock, and J. Eisert, Phys. Rev Lett. 100, 140409 (2008). 

[6] M. Cramer, Phys. Rev Lett. 106, 215302 (2011). 

[7] D.-S. Luhmann, K. Bongs, K. Sengstock, and D. Pfannkuche, 
Phys. Rev Lett. 101, 050402 (2008). 

[8] R. M. Lutchyn, S. Tewari, and S. Das Sarma, Phys. Rev A 79, 
011606 (2009). 

[9] A. Mering and M. Fleischhauer, Phys. Rev. A 83, 063630 
(2011). 

[10] A. Albus, R lUuminati, and J. Eisert, Phys. Rev A 68, 023606 
(2003). 

[1 1] C. Kollath et al, Phys. Rev A 69, 031601 (2004). 
[12] J. E. Hirsch, Physica C 158, 326 (1989). 
[13] R. Strack and D. Vollhardt, Phys. Rev Lett. 70, 2637 (1993). 
[14] J. E. Hirsch, Physica B 199/200, 366 (1994). 
[15] J. C. Amadon and J. E. Hirsch, Phys. Rev B 54, 6364 (1996). 
[16] G. Mazzarella, S. M. Giampaolo, and F. lUuminati, Phys. Rev. 
A 73, 013625 (2006). 



[17] D.-S. Luhmann, O. Jiirgensen, and K. Sengstock, New J. Phys. 

14, 033021 (2012). 
[18] U. Bissbort, F. Deuretzbacher, and W. Hofstetter, 

arxiv:1108.6047vl (2011). 
[19] J. Li, Y. Yu, A. M. Dudarev, and Q. Niu, New J. Phys. 8, 154 

(2006). 

[20] J. Larson, A. Collin, and J.-P. Martikainen, Phys. Rev. A 79, 
033603 (2009). 

[21] K. R. A. Hazzard and E. J. Mueller, Phys. Rev A 81, 031602 
(2010). 

[22] O. Dutta, A. Eckardt, R Hauke, B. Malomed, and M. Lewen- 

stein. New J. Phys. 13, 023019 (2011). 
[23] K. Sakmann, A. L Streltsov, O. E. Alon, and L. S. Cederbaum, 

Phys. Rev Lett. 103, 220601 (2009). 
[24] K. Sakmann, A. I. Streltsov, O. E. Alon, and L. S. Cederbaum, 

Phys. Rev. A 82, 013620 (2010). 
[25] T. Busch, B.-G. Englert, K. Rz^azewski, and M. Wilkens, Foun- 
dations of Physics 28, 549 (1998). 
[26] R R. Johnson, E. Tiesinga, J. V. Porto, and C. J. WiUiams, New 

J. Phys. 11, 093022 (2009). 
[27] S. Will, T. Best, U. Schneider, L. HackermuUer, D.-S. 

Luhmann, and I. Bloch, Nature (London) 465, 197 (2010). 
[28] H. R Buchler, Phys. Rev. Lett. 104, 090402 (2010). 
[29] S. Pilati and M. Troyer, arxiv: 1108. 1408 (2011). 
[30] S. Will, T. Best, S. Braun, U. Schneider, and L Bloch, Phys. 

Rev Lett. 106, 115305 (2011). 
[31] The resulting many -body product basis has a total length of 

12000n^, where n is the number of bosons. 



